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ABSTRACT 

Infrared stellar photometry from the Two Micron All Sky Survey (2MASS) and spec- 
tral line imaging observations of 12 CO and 13 CO J = 1-0 line emission from the Five 
College Radio Astronomy Observatory (FCRAO) 14m telescope are analysed to assess 
the variation of the CO abundance with physical conditions throughout the Orion A 
and Orion B molecular clouds. Three distinct A v regimes are identified in which the 
ratio between the 13 CO column density and visual extinction changes corresponding 
to the photon dominated envelope, the strongly self-shielded interior, and the cold, 
dense volumes of the clouds. Within the strongly self-shielded interior of the Orion A 
cloud, the 13 CO abundance varies by 100% with a peak value located near regions of 
enhanced star formation activity. The effect of CO depletion onto the ice mantles of 
dust grains is limited to regions with A v > 10 mag and gas temperatures less than 
~20 K as predicted by chemical models that consider thermal-evaporation to desorb 
molecules from grain surfaces. 

Values of the molecular mass of each cloud are independently derived from the 
distributions of A v and 13 CO column densities with a constant 13 CO-to-H2 abundance 
over various extinction ranges. Within the strongly self-shielded interior of the cloud 
(A v > 3 mag), 13 CO provides a reliable tracer of H 2 mass with the exception of the cold, 
dense volumes where depletion is important. However, owing to its reduced abundance, 
13 CO does not trace the H2 mass that resides in the extended cloud envelope, which 
comprises 40-50% of the molecular mass of each cloud. The implied CO luminosity to 
mass ratios, M/Lco, are 3.2 and 2.9 for Orion A and Orion B respectively, which are 
comparable to the value (2.9), derived from 7-ray observations of the Orion region. Our 
results emphasize the need to consider local conditions when applying CO observations 
to derive H2 column densities. 

Key words: ISM: abundances - ISM: clouds - ISM: individual - ISM: molecules - 
radio lines: ISM 



1 INTRODUCTION 

The mass and column density of interstellar molecular 
clouds impact their evolution, chemistry, and the produc- 
tion of newborn stars within their domain. Therefore, ac- 
curate measures of these key properties are essential to 
our understanding of this important interstellar gas phase. 
The millimeter and submillimeter rotational transitions of 
carbon monoxide ( 12 CO) and its isotopologues, 13 CO and 
C ls O, are the most widely used tracers of molecular hydro- 
gen (H2) that offer several advantages over other methods 
such as dust extinction, dust emission, and gamma rays. 
These include the ability to construct high angular resolu- 
tion views of distant interstellar molecular clouds through- 
out the Milky Way and other galaxies, and with sufficient 
spectral resolution, the provision of critical kinematic infor- 
mation such as the velocity centroid and velocity disper- 



sion for a given line of sight. Measurements of the opti- 
cally thick 12 CO J=l-0 line are routinely ma de in investiga- 
tions of star formati o n processes in galax ies (|Young et al. I 
Il995l : lOmont 1 120071 ; iLerov et alD l2009f) . The application 
of this high opacity line as a tracer of molecular hydro- 
gen relies on the empirical correlation between 12 CO ve- 
locity integrated intensity, Wco, and 7-ray intensity that 
is produced from the interaction of c osmic rays with hydro- 
gen protons within the Milky Way ( Bloemen et al. 1 1 19861 : 
IStrong fc Mattox 1 1 19961 : iHunter et al. 1 1 19971 ). The corre- 
sponding CO to H2 conversion f actor, Xco=N(H2)/Wc o 
is 1.9xl0 20 cm" 2 (K km s" 1 )" 1 (|Strong fc Mattox I(l996h . 
A similar procedure has been followed using the dust 
emission to trace gas column density with derived values 
of Xco ranging from 1.8-2.5xl0 20 cm" 2 (K km s' 1 ) ~ 1 
jDame. Hartmann. fc Thaddeus Il200ll ; lAde et al. Ikoill ) . 

Emission from 13 CO offers another measure of the 
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molecular gas distribution. Owing to its lower opacity, it 
is straightforward to derive the 13 CO column density under 
the assumption of local thermodynamic equilibrium (LTE), 
using the optically t hick 12 CO line to estimate the exci- 
tation temperature (jDickman I 1 19781 ). However, to derive 
the more relevant H2 column density and mass, knowl- 
edge of the 13 CO abundance is required. Most studies 
adopt a constant fractional abundance of ~2x 10 -6 through- 
out the cloud based on previous investigations of nearby 
dark clouds, which examined the relati onship between 
13 CO column density and visual ex tinction (jDickman II 19781 ; 
iFrerking. Langer. fc Wilson | [l982 V These earlier studies in- 
vestigated the N( 13 CO)-A„ relationship in cold, dark clouds 
with limited star formation activity. 

In this paper, we examine the relationship between car- 
bon monoxide and near infrared extinction in the Orion 
molecular cloud that is producing massive stars and presents 
a broad range of environmental conditions. In Section 2, we 
describe the near-infrared and molecular line data used in 
this study. The procedures used to derive extinctions from 
the 2MASS colours and the 13 CO column densities from the 
molecular line imaging are summarized in Section 3. In Sec- 
tion 4 the results of point-to-point comparison of the 13 CO 
column density and 12 CO integrated intensity to visual ex- 
tinction relation are presented. The interpretation of these 
relationships within the context of photon-dominated re- 
gions (PDRs) models, the calculation of relative abundances, 
and the derivation of cloud masses are described in Section 
5. 



2 DATA 

2.1 Millimeter Spectroscopic Imaging 

Observations of 12 CO and 13 CO J=l-0 emission were car- 
ried out with the FCRAO 14m telescop e using the 32 pixel 
SEQUOIA array (|Erickson et al. 1 1 19991 ) between 2005 and 
2006. Owing to the broad bandwidth of the HEMT am- 
plifiers, the 12 CO and 13 CO lines were observed simulta- 
neously enabling excellent positional registration and cali- 
bration. The backends were comprised of 64 autocorrelation 
spectrometers each with 25 MHz bandwidth. No smoothing 
was applied to the autocorrelation function so the spectral 
resolution was 29.5 kHz per channel corresponding to 0.077 
km s _1 and 0.080 km s _1 at the line frequencies of 12 CO 
(115.271202 GHz) and 13 CO (110.201353 GHz) respectively. 
System temperatures rang cd from 350-500 K ( 12 CO) and 
150-300 K ( 13 CO). The FWHM beam size of the antenna 
at the observed frequencies are 45" (115 GHz) and 47" (110 
GHz). The main beam efficiencies at these frequencies are 
0.45 and 0.48 respectively as determined from measurements 
of Jupiter. 

The areas of the Orion A and Orion B clouds were di- 
vided into submaps. A given submap of 12 CO and 13 CO 
emission was constructed by scanning the telescope at an 
angle 9 S = 14.04 degrees relative to the azimuth direction. 
This choice of scanning angle ensures even coverage of the 
sky along a single scan. Subsequent scans are orthogonally 
offset by the entire width of the array to extend the sky cov- 
erage. We also employed scanning angles of 6 S + 90. A basic 
observational block was comprised of a sequence of scans 



either using 6 a or 6 a + 90. Submaps were observed so that 
adjacent blocks overlap. The scanning angle for each block 
was set to be opposite to that of its neighboring block. The 
angular coverage spanned 22.2 deg 2 for Orion A and 17.4 
deg 2 for Orion B. The la errors due to thermal noise (main 
beam temperature units) within 0.2 km s _1 wide channels 
are 2.0 K (Orion A) and 2.1 K (Orion B) for 12 CO and 0.77 
K (Orion A) and 0.85 K (Orion B) for the 13 CO data. 

Images of velocity integrated 12 CO and 13 CO J=l-0 
emission from the Orion A and Orion B clouds at the full an- 
gular resolution are shown in Figure [T] and Figure [2] respec- 
tively. The irregularly shaped boundaries result from the 
scanning method. Data near the boundaries exhibit higher 
noise owing to less integration time. The images show sev- 
eral common emission features typical of giant molecular 
clouds. Regions of bright, high contrast emission are typi- 
cally sites of the most recent star formation. Such regions 
are often warm (T& > 25 K), filamentary, and character- 
ized by high column densities ( N(i?2) > 10 22 cm" 2 ). In 
addition, there is an extended, low surface brightness com- 
ponent that comprises most of the projected area of the 
cloud. Recent wide field imagery of 12 CO and 13 CO J=3- 
2 emi ssion from Orion A show s similar, large scale fea- 
tures (|Buckle et al. 1 120091 . |2012ft . Our wide field data on 
Orion B offers a significant improvement in sensitivity, an- 
gular arid_^rje£tralresohitiono^ previous maps of CO emis- 
sion l|Ballv. Langer. fc Liu lli"99ll ) 

2.2 Near-infrared Photometry 

Near-infrared extinctions are determined from colours of 
stars listed in the Point Sour ce Catalog of the Tw o Mi- 
cron All-Sky Survey (2MASS) (|Skrutskie et al. 1120061 ). The 
observations are comprised of three near-infrared bands: J 
(1.24/mi), H (1.66/im), and K s (2.16/im). The 2MASS Point 
Source Catalog achieves 99% differential completeness limits 
of 15.8, 15.1, and 14.3 magnitudes for J, H, and Ks bands re- 
spectively with 5% photometric accuracy for bright sources 
and 0"5 positional accuracy relative to the International Ce- 
lestial Reference System (ICRS). Positions and photomet- 
ric information are extracted from the 2MASS Point Source 
Catalog for objects that coincide with the area covered by 
the molecular line observations shown in Figure [T] and Fig- 
ure H The Orion A and Orion B fields contain 164,930 
and 242,931 sources respectively. In addition, six circularly 
shaped control fields with radius of 1° are analysed. These 
fields are selected based on their large angular displacement 
from known regions of molecular emission. The (I, b) centre 
positions of each control field are: (204, -22), (220,-20), (212,- 
23), (214,-16), (215,-23), and (220,-16). Only stars with de- 
tections in all three bands and with photometric errors less 
than 0.03 magnitudes are used for the control field, resulting 
in a total of 26,927 stars. 



3 ANALYSIS 

3.1 Gas Excitation and 13 CO Column Density 

Methods to derive column densities of 13 CO from ob- 
servations of the low rotational quantum numbers un- 
der the assumption of Local Thermodynamic Equilibrium 
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Figure 1. Images of (top) 12 CO and (bottom) 13 CO J=l-0 emission from the Orion A molecular cloud integrated over Vlsr to 16 
km/s at the full resolution of the observations. Units of the colour wedges are Kkms~ 1 . 



(LTE) are provided by iDickman I (|l978l ) and more recently, 
I Pineda et al. I j20TCljj. In brief, the excitation temperature 
that characterizes the population levels of the lower and up- 
per rotational states of 13 CO is determined from the peak 
brightness temperature of the optically thick 12 CO line, 
7i2,pfc- For the J=l-0 transition, 



5.53 



Ml + 



(1) 



Ti2,pfc +0.83 / 



The opacity of the 13 CO emission, 7"i3,i(«), in the upper, 
J=l rotational state is derived from the 13 CO brightness 
temperature, Tis(v), 
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1 There is a typo in equation 18 of Pineda ct al. (2010). The 
argument to the exponential function should read J (J + l)/kT x . 



The column density of 1 CO within the J=l state, Ni, is 
given by the expression 



iVi = 



8irkv 2 
hc 3 A 1{ 
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where Aio is the Einstein A coefficient and v is the frequency 
of the J=l-0 transition. To account for column density in the 
other rotational states, the assumption of LTE is applied, 

Z 
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where B is the rotational constant for 
partition function, 
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CO, and Z is the 
(5) 



The uncertainties in the 13 CO column densities are 
dominated by the assumptions of LTE. For regions with 
densities less than the critical density, the 13 CO rotational 



© 0000 RAS, MNRAS 000, 000-000 



4 F. Ripple, M.H. Heyer, R. Gutermuth, R.L. Snell, CM. Brunt 



120 



CD 
CD 
T3 





208 206 
I (deg.) 



204 



208 206 
I (deg.) 



204 



Figure 2. Images of (left) 12 CO and (right) 13 CO J=l-0 emission from the Orion B molecular cloud integrated over Vlsr to 16 km/s 
at the full resolution of the observations. Units of the colour wedges are Kkms -1 . 



energy levels may be subthermally excited. Thus, the excita- 
tion temperature derived from the more thermalized 12 CO 
J=l-0 line may overestimate the excitation temperature of 
13 CO. To illustrate this effect, we solve for the population in 
the 13 CO rotational energy levels using an escape probabil- 
ity radiative transfer method with the following conditions: 
volume density of 750 cm -3 , H2 column density of 3xl0 21 
cm~ 2 , a kinetic temperature of 10 K, a 13 CO to H2 abun- 
dance of 2xl0~ 6 , and a velocity dispersion of 0.4 km s _1 . 
The model excitation temperatures for the J=l-0 transition 
are 9.8 K for 12 CO and 5.5 K for 13 CO. In this case of sub- 
thermal excitation, the optical depth correction for 13 CO is 
underestimated leading to a smaller column density in the 
J=l level. This error is partially compensated by the over- 
correction for population of the higher excited states using a 
higher excitation temperature. For a spectrum with a peak 
antenna temperature of 1.5 K, the 13 CO LTE column den- 
sity derived using the 12 CO excitation temperature of 9.8 K 
is 0.7 times the column density using T ex —5 K. 

To assess the impact of the thermal noise inherent 
within the observations, we derive the excitation tempera- 
ture and 13 CO column density using a Monte Carlo method 
that estimates the statistical errors for each quantity. The 
calculations are performed on data cubes of the 12 CO and 
13 CO J=l-0 emission that are re-sampled to 1.8' pixel size 
and angular resolution. This resolution element corresponds 
to the pixel size of the infrared extinction map described 
in §3.2. For each location on the sky, values of peak 12 CO 
brightness temperatures and 13 CO integrated intensities are 
extracted from Gaussian distributions centered on the mea- 
sured value of the peak temperature of the 12 CO spectrum, 
Ti2,pfc, and W( 13 CO) = J Tn(v)dv with dispersions of a\2 
and CT13V AvSv where a 12 and an are the root mean-square 
values of antenna temperatures in spectral channels with 



no emission, Av is the channel spacing for each spectrum 
and Sv is the velocity width over which the 13 CO spectrum 
is integrated. In this study, we consider the Vlsr interval 
of to 16 km s _ so 5v=16 km s" 1 . A given set of Ti2, P k 
and W( 13 CO) values are propagated through equations 1- 
5 to produce realizations of T ex and N( 13 CO). This step 
is repeated 1024 times to generate distributions of T ex and 
N( 13 CO) values for each pixel. The value of the excitation 
temperature and statistical error at a given position, x,y, 
is respectively, the mean and dispersion of the T ex distri- 
bution. Similarly, the mean and dispersion of the N( 13 CO) 
distribution are assigned to the 13 CO column density and its 
statistical error respectively. Images of the excitation tem- 
perature and 13 CO column density are shown in Figure [3] 
for Orion A and Figure [1] for Orion B. Median 1-sigma un- 
certainties for T ex and N( 13 CO) in Orion A are 0.59 K and 
5x 10 14 cm -2 respectively. For Orion B, these values are 0.55 
cm' 



3.2 Visual Extinction 

The 2MASS Point Source Catalog is a powerful set of data 
to construct extinction images of nearby interstellar clouds. 
As an all sky survey, it covers the wide area subtended by 
nearby clouds. In addition, a greater number of background 
sources are detected in the near infrared relative to optical 
bands owing to the reduced opacity with increasing wave- 
length. To derive extinctions to each star, we implement 
the Near-Infra red Color Excess Revisite d (NICER) method 
described by lLombardi fc Alvesl (|2001). This method is a 
multi-band, maximum-likelihood approach for determining 
extinction to a background star based on the near-infrared 
colours. This is achieved through comparison of the observed 



stellar colour, 



with the intrinsic colour, c! n derived 
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Figure 3. Images of (top) gas excitation temperature, (middle) 13 CO column density, and (bottom) infrared derived visual extinction 
after applying the clipping algorithm for Orion A. Pixels for which the extinction is not defined are blanked and shown as white coloured 
pixels. These are located in the high column density regions of the cloud. The vertical dotted lines and numbers in the T ex image denote 
the partitions discussed in the text. 



from stars within a control field assumed to have no red- 
dening. With the three near infrared bands of 2MASS, two 
colours: ci = J — H and ci — H — Ks, are calculated. 
The inferred V-band extinction is derived from the mea- 
sured colour excess as 

(„obs _ int\ 

A Vti = lC ' - ' ' (6) 
hi 

where k, represents the ratio between the colour excess and 
the V band extinction. W e adopt the reddening law from 
Rieke & Lcbofskv (1985), corresponding to fci = 1/9.35 and 
k 2 = 1/15.87. The NICER method has the advantage of 
using any number of colours depending on the data available. 
Extinctions are derived in cases where photometric data are 
only available in two of the three bands. 

The intrinsic stellar J-H and H-Ks colours and colour 
variance of background stars are derived from stellar pho- 
tometry within the control fields. The mean colours and 



standard deviations for stars within the control fields are: 
0.446 ± 0.158 mag and 0.112 ± 0.058 mag for (J-H) and 
(H-K S ). 

For each star in the target field, the NICER procedure 
produces an extinction, A v , and its variance, Var(A„), that 
considers both the variance of colours in the control fields 
and the photometric errors of the target star. The uncer- 
tainty of these extinction estimates, (ta„ = y/Var(A v ), is 
dominated by the variance of the intrinsic colours of control 
field stars rather than the photometric errors of the target 
star. The median of derived extinction uncertainty values is 
1.4 magnitudes for both the Orion A and Orion B clouds. 

There are several sources of contamination to identify 
and remove prior to folding the derived extinction values 
into an image to compare with the 13 CO column density 
image. As the Orion cloud is an active region of star forma- 
tion, one expects a number of young stellar objects (YSOs) 
at various stages of pre-main sequence evolution within the 
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Figure 4. Images of (left) gas excitation temperature, (middle) 13 CO column density, and (right) visual extinction for Orion B. 



observed fields. Class I and II YSOs are associated with large 
infrared excesses due to the presence of a dense envelope 
or circumstellar disk that both attenuates the photospheric 
emission of the star and contributes to the near infrared flux 
from warm, thermally emitting dust grains. Therefore, the 
extinction inferred from the infrared excess of Class I and 
Class II YSOs is not representative of the column density 
through the ambient cloud. In Orion A, the degree of con- 
tamination by YSOs is particularly acute near the Orion 
Nebula C luster. The cata l og of Class I and II YSOs gen- 
erated bv lMegeath et al~l l|2012l ) derived from longer wave- 
length photometry from the Spitzer Space Telescope, is used 
to remove such objects from the list of background stars. A 
total of 2,456 (Orion A) and 553 (Orion B) Class I and Class 
II YSOs are excluded from our analysis. 

More problematic sources of contamination are fore- 
ground stars and Class III YSOs within the cloud. Both 
stellar types have intrinsic colours comparable to those de- 
rived in the control fields but whose light is not extincted 
by the full path of the cloud, if at all. The effect of includ- 
ing such stars in deriving the local extinction is illustrated 
in Figure [5] which shows a rare but extreme example of 
such contamination. A v values derived for stellar sources 
contained within a 1.8' xl.8' area of the cloud are rank- 
ordered for clarity. The low extinctions determined for 5 
stars is discrepant by ~ 20 magnitudes with the large A v 
values of four stars as well as the equivalent extinction based 
on the 13 CO column d ensity, A v>Bq =N( 13 CO)/2.0 x 10 15 
mag" 1 |Dickman Ill97&t ). While such extinction differences 
could be due to intrinsic, cloud structure within the area 
of the pixel, we attribute these discrepancies to contami- 
nation by foreground stars and Class III pre main sequence 
stars that are surely present within the Orion cloud that has 
been actively forming stars for several million years. Using 
extinction maps of star forming c louds at varying distances, 
lLombardi. Alves. fc LadlH (|201ll ) estimate 279 foreground 
stars per square degree towards the Orion cloud. This corre- 
sponds to 0.25 foreground stars per 1.8' sized pixel. Owing 
to magnetic activity in their atmospheres, Class III YSOs 
can be identified by their X ray emission. A list of X-ray 
identified Class III stars in Orion A (L1641 and NGC 1980) 



has b een compiled from XMM-Newton data dPillitteri et al. I 
l2012h . In total, 443 Class III stars from this list have 2MASS 
counterparts and are excluded from our data analysis. 

The X ray survey of the Orion A cloud does not pro- 
vide a complete census of Class III YSOs nor does it iden- 
tify foreground stars. To isolate and remove additional Class 
III YSOs and foreground stars, we apply the following al- 
gorithm that conservatively leverages the additional infor- 
mation provided by the 13 CO column density. If A v , eq is 
greater than an extinction threshold, A V)t h, and the extinc- 
tion of any star within the 1.8' pixel is less than a frac- 
tion of the equivalent extinction, fth.Av,eg, then the star is 
excluded from any subsequent analysis. Optimal values for 
A v .th and fth are determined by maximizing the number 
of X-ray emitting Class III YSOs identified by the algo- 
rithm. From this optimization, we find K Vi th = 8 mag and 
fth = 30%, which identifies 56% of the 443 known Class III 
stars. The algorithm is illustrated in Figure^ The solid line 
shows the A v , eq estimated from N( CO) and the dotted 
line shows A v ,th below which our algorithm would remove 
stars. Following the removal of the contaminating stars, the 
mean extinction shifts from 8.7 (dashed line) to 18.2 magni- 
tude (dotted-dashed line). A more liberal cut by increasing 
A v .th and/or f t h would increase the likelihood of removal of 
valid background stars. The algorithm is limited to regions 
with large 13 CO column densities. Applying the algorithm 
to all stars, 2090 additional, potentially contaminating stars 
(1.1% of total) are isolated and removed from the list of 
background stars for Orion A. For Orion B, 824 contami- 
nated sources are culled from this list. Contamination by 
foreground and /or embedded sources may still be present in 
regions of lower 13 CO column density. 

An image of dust extinction is created from the set of 
infrared colours of background stars within our final list fol- 
lowing the removal of YSOs and foreground stars. Spatial 
grids are setup that are aligned with the respective images 
of 13 CO column density with 1.8' size pixels. The extinction 
of all sources whose position falls within the solid angle sub- 
tended by the pixel are averaged with equal weighting. The 
resultant image has a fixed resolution throughout the map 
at the expense of variable noise owing to a varying number 
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Figure 5. The rank-ordered values of derived extinction for 9 stars contained within a 1.8x1.8' area. The effect of contaminating 
foreground and diskless YSO sources is to lower the mean extinction. The heavy solid line represent the A v value determined from 13 CO, 
the dotted line shows the threshold, A v t ^ of the clipping algorithm, the dashed lined denotes the mean extinction value before the cut, 
and the dotted-dashed line shows the mean extinction value after removing the contaminating sources. 



of extinction values that are averaged. The median la ex- 
tinction errors for these sized pixels for A v < 10 are 0.8 mag 
and 0.7 mag for Orion A and Orion B respectively. For A v 
^ 10, the median extinction error is 1 mag for both clouds. 
These values are the propagation of the random errors for 
the set of extinctions that are averaged into a pixel. The 
standard deviation of the set of extinction values is typically 
larger than the associated random error. For extinctions less 
than 10 mag., the median dispersions are 1.2 and 1.4 magni- 
tudes for Orion A and Orion B respectively. For A v > 10, the 
dispersions are 2.5 and 2.3 magnitudes. In high extinction 
regions, there are many pixels through which no background 
stars shine. In these cases, the pixel is blanked and excluded 
from all subsequent analyses. The number of blanked pixels 
in Orion A and Orion B are 1475 and 1002 respectively. The 
final extinction maps constructed from the set of infrared 
extinctions are shown in Figure [3] and Figure [4] The me- 
dian extinction values over the entire map for Orion A and 
Orion B are 2.1 and 2.2 magnitudes respectively. For both 
clouds, the extinction map resembles the basic outline of the 
13 CO column density maps. 



4 RESULTS 

4.1 N( 13 CO)A„ Relationship 

Our analysis of near infrared stellar photometry from 
2MASS and FCRAO molecular line data generates two inde- 



pendent views of the mass distribution in the Orion molec- 
ular cloud. Dust extinction provides an excellent tracer of 
hydrogen column density, N(H) = N(HI) + 2N{Hi), for 
a given line of sight based on the established relationships 
from UV spectroscopy, N(H) = 5^xlCF E(B - V) cm" 2 



mag^ 1 (jBohlin. Savage, fc Drake"! Il978l: 

1 (|V 



|200Sh and A„/E(B — V) = 3.1 (|Whittet 



lachford et al. I 
20031 ). The "hT 



column density is derived from the 13 CO column den- 
sity with knowledge of the 13 CO abundance, N(-Ha) = 
N( 13 CO)/[ 13 CO/H 2 ]. However, the 13 CO abundance is ex- 
pected to vary with the local environment conditions due to 
changing UV radiation field intensity, gas density structure, 
and chemical depletion. 

Using the extinction as a reliable measure of the H2 col- 
umn density, we assess the 13 CO abundance by examining 
the relationship between N( 13 CO) and A v . The variations 
of 13 CO column density with A v for the entire Orion A and 
Orion B clouds are displayed in Figure [6] as shaded 2 di- 
mensional histograms within bins of A v and N( 13 CO). For 
the zoomed figures, the solid circles denote median values 
of N( 13 CO) derived within varying sized bins of A v which 
contain 250 points. The vertical and horizontal error bars 
reflect the standard deviation of N( 13 CO) and A v values re- 
spectively in each bin. Different trends are evident within 
three distinct regimes of extinction. Within the first regime 
at low A v (0 < A v < 3 mag), N( 13 CO) grows very slowly 
with increasing A v . The second extinction regime (3 ^ A v 
^ 10 mag) is characterized by much steeper linear growth of 
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Figure 6. (left) 2D histograms showing the variation of N( 13 CO) with A v for all defined pixels in the Orion A (top) and Orion B 
(bottom) clouds over the full range of visual extinction, (right) An enlarged version of the same figure over a more limited A v range. The 
solid circles and vertical error bars represent the mean N( 13 CO) values and standard deviation respectively within equally populated A v 
bins. Horizontal error bars show the dispersion of A v values within each bin. The white and black solid lines show the fit of equation 7 
for A v < A v fi and equation 8 for j4„ q < A v < 10 respectively. 



N( 13 CO) with A v . The third regime at high extinctions (A v 
> 10 mag) shows a continued linear trend of N( 13 CO) but 
with a much larger scatter. For Orion A, the large scatter is 
due, in part, to a secondary set of points in which the 13 CO 
column density remains constant at the value ~ 2.5xl0 16 
cm -2 with increasing A v . 

To quantify this transition from the low to interme- 
diate A v regime, a non-linear least squares (Levenberg- 
Marquardt) algorithm is applied to simultaneously fit two 
slopes (mi and m^) and an intersection point (A Vt i n t,Ni n t) 
to the data over the range < A v ^ 10 magnitudes, 

N{ 13 CO) = mi{A v - A Viint ) + N lnt (7) 

for < A v < A v ,i„t and 

N( 13 CO) = m 2 (A v - A v ,int) + N mt (8) 

for A v ,i nt ^ A v < 10. This piecewise linear equation takes 
into account the significantly different N( 13 CO)-.A„ relation- 
ship of the first two regimes. The fitted parameters of equa- 



tions 7 and 8 for Orion A and Orion B are listed in Table Q] 
and the fits are shown in Figure [6] These values describe the 
observed trends in Figure [6] in which N( 13 CO)/A„ ratio is 
small for A v < A v ,i n t and increases by factors of ~5 (Orion 
A) and ~8 (Orion B) for A v > A v ,mt- 

The Orion A cloud offers a wide range of physical con- 
ditions (FUV radiation field, gas temperature, column den- 
sity) that can impact the relation between N( 13 CO) and 
Av To assess the effect of the local environment on the 
N( 13 CO)-^„ relation, we divide the cloud into 9 Galactic 
longitude partitions, each containing ~2,500 pixels of non- 
blanked data (see Figure [3j . The Orion Nebula, M43, and 
NGC 1977 are located in partition 6 and make up part of 
the larger Orion OBI Association that extends from parti- 
tion 5 through partition 9. The N( 13 CO)-A„ relation for each 
partition is shown in Figures 09] In these cases, the solid 
circles and error bars in the zoomed figures reflect the me- 
dian and standard deviation respectively of N( 13 CO) values 
in A v bins containing 50 points. The trends in the three A v 
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Table X. N( 13 CO) - A v Relationship Fitted Parameters 



Partition mi m 2 A-v,int N int A„ j3 

(10 15 cm" 2 mag -1 ) (10 15 cm -2 mag -1 ) (mag) (10 15 cm -2 ) (mag.) 



Orion B 


0.42 (0.05) 


3.19 


(0.20) 


2.98 (0.05) 


1.44 (0.07) 


3.5 (0.1) 


Orion A 


0.47 (0.04) 


2.17 


(0.08) 


2.68 (0.04) 


1.84 (0.06) 


3.2 (0.1) 


1 


0.4 (0.2) 


1.5 


(0.1) 


2.7 (0.2) 


3.0 (0.2) 


2.7 (0.2) 


2 


0.7 (0.1) 


2.3 


(o.i) 


2.7 (0.1) 


2.7 (0.2) 


2.8 (0.1) 


3 


0.4 (0.2) 


1.8 


(0.1) 


2.2 (0.1) 


2.0 (0.1) 


2.8 (0.1) 


1 


0.5 (0.1) 


4.4 


(0.3) 


3.2 (0.1) 


2.4 (0.1) 


3.3 (0.1) 


5 


0.8 (0.3) 


2.0 


(0.2) 


2.3 (0.2) 


2.2 (0.3) 


2.7 (0.3) 


6 


0.5 (0.1) 


3.6 


(0.3) 


3.7 (0.1) 


2.1 (0.2) 


4.0 (0.1) 


7 


0.3 (0.1) 


1.7 


(0.4) 


3.1 (0.2) 


1.2 (0.1) 


4.2 (0.3) 


8 


0.1 (0.1) 


0.9 


(1.2) 


2.9 (0.8) 


0.9 (0.2) 


5.2 (3.2) 


9 


-0.07 (0.1) 


0.6 


(0.5) 


2.4 (0.3) 


0.7 (0.1) 


6.2 (3.2) 
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Figure 7. Same as Figure \E\ for partitions 1 (top), 2 (middle), and 3 (bottom) of the Orion A cloud. 
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regimes observed in Figure [6] are evident in these localized 
regions of the cloud. For each partition, we fit Equations 7 
and 8 to the 50 point A v bins for A v ^i 10 magnitudes. The 
results from these fits are summarized in Table [T] 

The fitted slopes in the low A v regime, mi, are similar 
for partitions 1 through 6, ranging from 0.4-0.8xl0 15 cm~ 2 
mag -1 . The northernmost partitions (7-9) show the lowest 
mi values, < 0.3x 10 15 cm -2 mag -1 . These values are con- 
sistent with little or no detected 13 CO in these areas. In the 
intermediate A v regime, the slope, m,2, increases by a fac- 
tor of 3 from partition 1 to partition 4. For partitions 4-6, 
the 2-line fit does not accurately describe the set of points 
in the range 5 < A v < 10 mag in which N( 13 CO) values 
are in excess with respect to the fitted line. In these cases, 
the fitted parameters are leveraged by the large number of 
points between A v ^ nt and 5 magnitudes with a lower ratio 
of N( 13 CO) and A v . For regions closer to the central region 
of the OBI association (partitions 7-9), the slopes rapidly 
decrease. As noted in Section 3.1, the measured 13 CO col- 
umn densities may be underestimated in the low volume 
density and low column density regions of the cloud owing 
to subthermal excitation of the J=l-0 transition. However, 
this error does not fully account for the measured differences 
between ml and m2 over the expected conditions within the 
cloud responsible for detectable 13 CO emission. 

The third A v regime also shows differential behavior. In 
partitions 2-5, N( 13 CO) does not continue to increase with 
increasing values of A v above 12 mag. of visual extinction. 
However, in partition 6, N( 13 CO) shows continued linear 
growth to visual extinctions of 30 mag. The different behav- 
ior in this partition accounts for the bifurcation of points in 
Figure [S] at high A v values. Partitions 1,7,8, and 9 lack suf- 
ficiently high extinction to define the behavior of N( 13 CO) 
in this regime. 



4.2 The Wco-A v Relationship 

Velocity-integrated 12 CO J=l-0 emission, Wco, is the most 
frequently applied measure of H2 column density and mass 
in galaxies owing to its brightness and correlation with 
other tracers of H2 such as 7 ray and dust emiss i on in 
the Milky Way jXde et al. I l201ll ; iBloemen et al. I 1 19861 ; 
IStrong fc MattoxlTl996h . Here, we assess the point-to-point 
relationship of Wco with the infrared derived extinction in 
the Orion A and B clouds and within the Orion A parti- 
tions. Figure [10] displays the variation of Wco with A v over 
the surveyed areas of the Orion A and B clouds. As a conve- 
nient reference, the line corresponding to the standard con- 
version of Wco to H2 column density, Nh2 = 1-9 x 10 20 Wco 
cm~ 2 /(Kkms~ 1 ) and N H 2 = 9.4 x lO 20 ^ cm -2 is also 
drawn on the plots. As found in previous studies, Wco 
correlates well with A v over the limited extinction range 
of ~2-5 magnitude s jPineda. Caselli. fc Goodman I 120081 ; 
iPineda et aTl l2010h . Yet for a fixed value of A v , there is 
a large scatter of Wco values. This scatter limits the use 
of 12 CO emission as an accurate measure of the H2 column 
density for a given line of sight through a molecular cloud. 
For extinctions less than 2 magnitudes, there is a large num- 
ber of points below the reference line corresponding to weak 
or no detected signal of 12 CO emission. These points corre- 
spond to gas layers in which there is significant dust extinc- 
tion but faint, underluminous 12 CO emission with respect 



to typical ratios of Wco to A v . Since the derived A v values 
are measured with respect to the extinction in the reference 
fields of background stars, the dust coupled to the atomic 
gas envelope has been largely removed. Most, if not all, of 
the extinction in the range of 0.1 to 2 magnitudes, is linked 
to dust associated with molecular hydrogen. 

To better understand the origin of these various trends 
in the Orion clouds, we examine the Wco~A v relation- 
ship within different environmental conditions. Figure [TTJ 
shows this relationship within each partition of the Orion 
A cloud. In the lower partitions (1-5), Wco linearly cor- 
relates with A v between the range 1-5 magnitudes. For 
partitions 2 and 3, Wco appears t o saturate for A v > 5. 
IPineda. Caselli. fc Goodman I (|2008l ') find a similar trend in 
the subregions of the Perseus molecular cloud and attribute 
the flattening of Wco with increasing A v to high optical 
depths rather than depletion of 12 CO onto dust grains. In 
the areas closer to the OB association (partitions 6-9), there 
are many positions with A v between and 4 magnitudes 
for which CO emission is anomalously faint. These points 
correspond to underluminous locations noted in Figure 10 
but can now be isolated within the cloud in the vicinity of 
the Orion OB stars. 



5 DISCUSSION 

In §4, we identify 3 distinct regimes in the relationship be- 
tween the 13 CO column density and the infrared-derived vi- 
sual extinction. In the following sections, we assume that the 
measured trends reflect the chemistry of 13 CO exposed to 
varying degrees of UV radiation with depth into the cloud, 
which in turn, depends on the local physical conditions of the 
cloud. The abundance of 13 CO relative to H2 is calculated 
using the derived 13 CO column densities and dust extinc- 
tion as a proxy of H2 column density. These abundances are 
examined in the context of photochemistry models and the 
changing physical conditions within the Orion cloud. Line 
excitation and optical depth of the 13 CO J=l-0 line can also 
impact the observed N( CO)-^,, relationships. In regions of 
low volume density, (n < 10 3 cm -3 ), that may characterize 
the low column density envelopes of GMCs, the J=l-0 tran- 
sition is subthermally excited, (T ex < kinetic temperature) 
leading to weak emission that is difficult to detect in low sen- 
sitivity, wide field mapping programs such as the molecular 
line data presented in this study. In regions of high column 
density, (A v > 15 mag.), the 13 CO J=l-0 emission may be 
optically thick and our method may not fully account for this 
effect. Our limited 12 CO and 13 CO J=l-0 data does not dis- 
tinguish such excitation and opacity effects from chemical 
variations. 



5.1 Baseline PDR Model 

To illustrate our observational results, the Meudon PDR 
model is used to compute the variation of 13 CO abundances 
under varying physic al conditions and UV radiation fields 
|Le Petit et al. Il2006h . The model describes the physical and 
chemical conditions within a slab of gas and dust with hy- 
drogen density, nn, and an incident FUV radiation field 
strength, x> measured in units of the iDraine (1978) radi- 
ation field. Although the model does not consider the de- 
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Figure 8. Same as Figure [6] for partitions 4 (top), 5 (middle), and 6 (bottom). 



pletion of molecules onto grain surfaces at high densities, 
the photochemical calculation is sufficient for comparisons 
to data in the low to intermediate A v regimes. The PDR 
models are calculated with parameters x = lj 10, 100, 10 3 , 
and 10 4 at a density of 10 3 cm -3 and x = 10 3 and 10 4 at a 
density of 10 4 cm -3 . These values of x an d volume density 
bracket the expected range of physical conditions through- 
out most of the Orion cloud as traced by 13 CO J=l-0 emis- 
sion. For all model parameter sets, the UV radiation field is 
assumed to be incident on both sides of the slab. It is impor- 
tant to note when comparing the model results to the ob- 
servations, that the model A v values are measures of depth 
into the cloud with the corresponding 13 CO column density 
to each depth whereas the observed A v and N( 13 CO) values 
shown in Figures [5121 reflect the total extinction and column 
density respectively through the cloud. For each depth po- 
sition in the slab as measured by A v , the model calculates 



the fractional abundance of atomic and molecular hydrogen 
and selected molecules with respect to the total density of 
hydrogen protons, nu = tihi + 2uh2- 

The N( 13 CO)-A„ relationships produced by the models 
with riH = 10 3 cm -3 are shown in Figure [12] These plots ex- 
hibit similar profiles as seen in Figures [BJH for the low and 
intermediate A v ranges — very small levels of 13 CO col- 
umn density within the low A v regime and a larger, yet con- 
stant value of N( 13 CO)/j4„ for the intermediate A v regime. 
For a given volume density, the depth or equivalent, model 
A v value above which the 13 CO abundance is constant, in- 
creases with increasing UV radiation field intensity (increas- 
ing x)- A higher radiation field strength requires a larger ex- 
tinction to fully self-shield 13 CO molecules to achieve a max- 
imum abundance. For cloud models with n=10 3 cm -3 , the 
depths at which 95% of the hydrogen is molecular are: A v — 
0.003, 0.12, 0.9, and 3.1 magnitudes for x=l, 10, 100,1000 re- 
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Figure 9. Same as Figure [6] for partitions 7 (top), 8 (middle), and 9 (bottom). 



spectively. At these depths the CO and CO abundances 
are several orders of magnitude smaller than their nomi- 
nal, strongly self-shielded values of lxlO -4 and 2xl0 -6 re- 
spectively. A consequence of the differential requirements for 
self-shielding of H2, 12 CO and 13 CO is a layer of molecular 
hydrogen gas within the cloud envelope that is not traced by 
CO owing to its negligible abundance. This cloud component 
of mol ecular hydrogen with no CO has been labelled "dark" 
gas bv lGrenier et al.l (|2005h and may comprise a significant 
mass reservoir that must be considered when estimating the 
self-gravity and surface density of GMCs. 



5.2 CO Abundance Variations from 

UV-Illuminated Cloud Envelope to the 
Self-Shielded Interior 



The low extinction regime identified in our analysis rep- 
resents areas within the cloud where the transition to 
fully self-shielded 13 CO material has not yet occurred such 
that photodissociation still dominates the chemistry of 
13 CO. The shallow, non-zero slope, mi, in this regime re- 
flects the low effective abundance of 13 CO. To derive the 
mean relative abundance from these data, [ 13 CO/H2] = 
<N( 13 CO)/yl l ,/(N(//2)/A„> = mi/9.4 x 10 20 . The values 
range from (1.1-8.5) x 10 -7 along the length of the Orion A 
cloud. The PDR models show that 13 CO abundance changes 
rapidly in this regime so these abundance estimates may not 
be representative of more narrow A v intervals. The models 
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Figure 10. Two dimensional histograms of A v and Wco values for (left) Orion A and (right) Orion B. The solid circles and error bars 
reflect the mean and standard deviation within A v bins containing 2048 individual points. The dashed horizontal line denotes the 3<r 
threshold of the Wco- The white line shows the relationship implied by the standard conversion between Wco an d H2 derived from 
7-rays (Strong et al. 1996). For both clouds, Wco i s underluminous with respect to this standard value for A v < 2 magnitudes. 



show that there is an appreciable amount of H2 mass within 
this layer. 

At sufficient depths into the cloud, self-shielding enables 
the formation of appreciable column density of CO. In this 
self-shielded regime, a steeper slope, 7712 of the N( 13 CO)-A„ 
relationship is measured, relative to slopes derived in the low 
A v regime and reflects the increased abundance of 13 CO. 
Within this regime, the relative abundance is [ 13 CO/H2] = 
rri2/9.4 x 10 20 . There are significant variations of the 13 CO 
abundance throughout the Orion A cloud. At the high lon- 
gitude end of the cloud (partitions 1-3), the average abun- 
dance is 2.0 x 10 -6 . In the central partitions of the cloud 
(partitions 4-6), the mean abundance increases to ~3.5x 
1(T 6 . From this peak, the 13 CO abundances decrease in par- 
titions 7 and 8. Finally, the measured slope in partition 9 is 
consistent with no appreciable 13 CO abundance. 

These abundance results for the Orion A and Orion 
B molecular clouds are compared to values derived from 
other interstellar clouds with different conditions. Using 
V band st ar counts t o esti mate A v towards 38 interstel- 
lar clouds, IDickmanl l| 19781 ) found a mean relative abun- 
dance of 2.0 x 10 -6 with a scatter of 50% over the range 
1.5 <A V < 5 magnitudes. With the availability of near in- 
frared data to measure extinction in high c olumn den- 
sity regions, iFrerking. Langer. fc Wilson I (| 19821 ) determined 
the 13 CO abundances of 1.5 xlO" 6 ( 1 <A V < 5) and 2.9 
xl0~ 6 (4 <A V < 15) for the Tau rus and p-Oph m olec- 
ular clouds respectively. Fo r both IDickmanl (|l978l ) and 
IFrerking. Langer. fc Wilson I (|l982T ). the abundance values 
are based on a limited number of positions in each cloud. 
The advent of large-scale, multiband near infrared observa- 
tions enabled more spatially complete est imates of mo l ecular 
abundances to deeper extinction regimes. lLada et al. I (| 19941 ) 
investigated the dark cloud IC 5146 to A v depths of 32 mag 
of extinction to derive a 13 CO abundance of 2.3 x 10 -6 . More 
recent studies of the N( 13 CO)-j4„ relationship have analysed 
2MASS data. iPineda. Caselli. fc Goodman I |2008l ) derived 
a mean 13 CO abundance of 2.5x 10~ 6 for the Perseus molec- 



ular cloud. Their study also found ~25% variations of the 
13 CO abundance with the largest ab undance values in the 
active star forming region, NGC 1333. IPineda et all (|20ld ) 
found a m ean abundance of 1.6xl0~ 6 in the Tau rus molec- 
ular cloud. lHariunpaa. Lehtin. fc Haikala 1 |2004) examined 
a diverse set of molecular clouds to find a broad range of 
13 CO abundances (0.9-3.5 x 10~ 6 ) with the largest values 
linked to active star forming regions. 

The conditions within the Orion A cloud span those 
found in both dark clouds as well as active star forming re- 
gions. The effective, global abundance of the Orion A cloud 
is similar to that of p Oph, which is a moderately active 
region of star formation containing a limited number of O 
and B stars. Partitions 1-3 are typical of quiescent star form- 
ing regions such as Taurus and IC 5146 and exhibit abun- 
dances similar to values found in these dark clouds. The 
abundances in the warmer, star-forming regions (partitions 
4-6) of Orion A are more comparable to those found in ac- 
tive, star-forming regions. T he partitions in Orion A show a 
simila r trend as described in lHariunpaa. Lehtin. fc Haikalal 
(2004). The much lower abundances found in partitions 7 
and 8 are likely a result of reduced extinction and a more 
intense UV radiation field. For partition 9, the second and 
third A v regimes are non-existent, implying that the UV 
field intensity is too large for 13 CO to become self-shielded. 

The elevated 13 CO abundance values in regions of ac- 
tive star forming regions within Orion and other targets 
pose an important constraint to descriptions of interstel- 
lar chemistry. Such enhanced abundances may result from 
the destruction of ice mantles on dust g rains and t he re - 
lease of previously depleted CO molecules. I Williams I (|l985h 
summarizes several mechanisms including the evaporation 
of the mantles from the energy released in radical reactions 
on the dust grains, photodesorption by UV radiation trig- 
gered by cosmic ray interactions with H2 molecules in the 
interior of the cloud, and shocks generated from the inter- 
action between dense clumps and stellar winds from young 
stellar objects. 
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Figure 11. Scatter plots of the A v and Wco values for each Orion A partition. The line shows the relationship implied by the standard 
conversion between Wco an d H2 and the dashed horizontal line denotes the 3cr threshold of the Wco- The overdensity of points for 
A v <3 and Wco < 4 K km s — 1 in partitions 6-9 reflect underluminous 12 CO emission owing to reduced 12 CO abundance in this regime. 



The transition of CO abundances from the UV- 
illuminated envelope to the self-shielded interior are depen- 
dent on the physical conditions. To quantify this transition 
and enable a comparison of the data with models, we evalu- 
ate the Av depth, A v ,3, at which N( 13 CO) equals a nominal 
yet significant value of 3xl0 15 cm -2 . For each partition in 
the Orion A cloud, A Vt z is calculated from the parameters, 
rri2, A v ,i nt , and Nint- Figure [13] shows the variation of this 



extinction threshold and model-derived thresholds for dif- 
ferent values of x/ n - The extinction thresholds are not well 
characterized in partitions 8 and 9 owing to the absence 
of any sig nificant 13 CO emission and therefore, are not in- 
cluded in this plot. A v ,3 is nearly constant at ~2.7 mag. in 
partitions 1-3 and partition 5 and is well bounded by models 
with 0.001 < x/ n <0.005. Given the large displacement of 
partitions 1-3 from the bright ionizing sources of the Orion 
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Figure 12. The relation between N( 13 CO) and A v derived from models produced by the Meudon PDR code for n=10 3 cm 3 and 
different values of the radiation field, y> (see legend). The model does not include the effects of 13 CO depletion onto dust grains. 



OB association, we estimate the conditions x ~ 1-5 arm n ~ 
10 3 cm -3 in these areas. Partition 5 lies closer to the OB 
association so one might expect an enhanced UV field but 
this could be balanced by a corresponding increase in the 
gas density. For partitions 6 and 7, A Vi 3 increases to 4 mag- 
nitudes and is well matched with model values with x/ n of 
~0.1. This condition can be reasonably matched to x=10 3 
and n=10 4 cm -3 . The absence of significant 13 CO emission 
in partitions 8 and 9, which lie near the centre of the OB as- 
sociation, likely indicates a very intense ultraviolet radiation 
field and a low density medium for which the condition for 
13 CO self-shielding is not possible. The measured changes 
in the threshold extinctions in the various partitions reflect 
the varying intensity of the far ultraviolet radiation field 
owing to displacement from the Orion OBI association and 
Orion Nebula Cluster and the local volume density condi- 
tions within the cloud. 



5.3 High A v Regime 

The role of depletion of molecules onto interstellar dust 
grains in cold, high density regions of molecular clouds has 
been increasingly recognized owing to the availability of 
secondary tracers of H2 column density that probe deep 
into the cloud, such as infrared derived extinct ion and sub- 
millim eter thermal dust continuum emission. lLada et al. I 
l|l994h observed a flattening of 13 CO integrated intensity 



in the I C 5146 cloud for lines of sight with A v > 8-10 mag- 
nitudes. iKramer et al. I fl999) found a similar trend in IC 
5146 for C ls O column densities for A v > 10, which they 
attributed to depletion . Other examples o f depletion in- 
clude the globule B 68 (|Bergin et al. I |2 002f ) , starless cores 
jTafalla et al. 1 12002| 1. filaments within the Taurus cloud 



l|Pineda et al. 



l2010fl . IChristie et all 1 1201 2) examined de- 



pletion effects in a sample dust cores within several Gould's 
Belt clouds, including regions in Orion B. They found de- 
pletion of C ls O to be most severe in the Taurus cloud and 
smallest in the Ophiuchus cloud. Both prestellar and pro- 
tostellar cores within the Orion B cloud exhibit varying de- 
grees of C 18 depletion. 

Models of molecular cloud chemistry that consider de- 
pletion ident ify several key factors that modulate molecular 
abundances. iBergin. Langer. fc Goldsmith I (1995) show de- 
pletion of CO depends strongly on dust temperature, den- 
sity, and timescale. At low densities (<10 4 cm -3 ) and ex- 
tinctions, the effect of depletion is small owing to desorption 
by cosmic rays and radiation that maintain significant frac- 
tions (60-90%) of CO in the gas phase, even at late times of 
10 7 years. For densities in excess of 10 5 cm -3 , dust tempera- 
tures of 10 K, and A v = 10 mag., the depletion time scales are 
short and much of the CO (~96%) freezes onto the grain sur- 
faces. For increasing grain temperature, molecules attached 
to grain surfac es are released by thermal evaporation of the 
grain mantles. B ergin. Langer. fc Goldsmith I jl995) show 
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Figure 13. Values of the extinction, A v .s, at which N( 13 CO)= 3xl0 15 cm -2 for partitions 1-7 in the Orion A cloud. The horizontal 
lines correspond to A v ^ calculated for each PDR model described by different values of xl n - The variation of A v ^) in the cloud is linked 
to proximity to the Orion OB association and changes in the local volume density of the gas. 



that thermal evaporation is an extremely effective desorp- 
tion process at dust temperatures greater than 22 K with the 
degree of chemical depletion depending on the effectiveness 
of the desorption processes. For temperatures larger than 
~25 K the models calculate that all CO is expected to be 
found in the gas phase However, temperature programmed 
desorption (TPD) laboratory studies demonstrate that CO 
molecules can be locked into water ice on grain mantles and 
that CO remains partially (25-40 %) depleted on these sur- 
faces to even higher temperatures l|Collhigs et al. |[2004l ). 

In the Orion A cloud, differential behavior in the high 
A v regime is identified. In partitions 2-5, the extinction at 
which the N( 13 CO)-A„ relation flattens represents the depth 
at which environmental conditions are appropriate for 13 CO 
depletion onto dust grains. The turn-over occurs at approx- 
imately 10 mag in partitions 2-4 and 15 mag in partition 5. 
In contrast, partition 6 shows continued linear growth with 
increasing A v above 10 mag. Owing to the limited range of 
A v in partitions 1,7,8, and 9, it is not possible to assess the 
role of depletion in these areas. 

To examine the role of temperature in the depletion pro- 
cess, the gas temperatures derived from the optically thick 
CO emission in each partition are examined as a proxy or 



relative measure to the grain temperature, considered in the 
models. Table [5] shows the mean excitation temperatures 
T ex for two regimes of extinction, moderate (5 < A v <10) 
and high (A v > 10), chosen to ensure that the CO lines are 
thermalized and to sample the high column density regions 
where depletion occurs. The temperature within the high 
extinction regions of partitions 2-4 are all below the 22 K 
threshold calculated in the models, while partition 5 reaches 
25 K. In partition 6, where the temperatures are significantly 
higher (T ex = 34.0 K), a significant fraction of 13 CO remains 
in the gas phase to A v =25 magnitudes. The limited amount 
of depletion in Partition 6 may also explain why it has the 
largest CO abundance in the cloud. Our data can not assess 
whether a significant amount of CO remains entrapped in 
water ice as measured in the TPD studies. If a large amount 
(50%) of CO remains depleted in these regions, as measured 
in the TPD experiments, then the interstellar abundance of 
CO is much larger than conventionally assumed. 

Thus, CO depletion is observed in the cold, dense re- 
gions of the Orion A cloud. The scatter and bifurcation of 
the Orion A composite N( 13 CO)-A„ relation shown in Fig- 
ure [6] is attributed to varying conditions within the cloud 
that allow chemical depletion of 13 CO onto dust grains in 
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Table 2. Properties within Partitions 
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the cold, dense regions and desorb 13 CO molecules in the 
warmer regions of the cloud. 



5.4 Cloud Masses 

An important goal of this study is to assess the reliability 
of 12 CO and 13 CO J=l-0 emission as a tracer of Eh mass 
within interstellar clouds under varying physical conditions. 
This effort is critical to studies of the molecular ISM and 
star formation throughout the Milky Way and within nearby 
galaxies, where infrared-derived extinctions over the angu- 
lar extent of GMCs are not possible. Previous calibrations 
of the 13 CO abundance in cold, dark clouds may not apply 
to the more massive GMCs that contribute the bulk of the 
CO luminosity in galaxies. The global masses of the Orion 
A and B clouds and the masses within the individual par- 
titions of Orion A are calculated for different cloud layers 
using the derived A v distributions, Hi 21 cm line emission, 
12 CO and 13 CO J=l-0 emission. For the atomic gas com- 
ponent, we analyse the 21 cm line emission from the Atla s 
of Galactic Neutral Hydrogen (|Hartmann fc Burton |[l997f ) . 
For all mass calculations, we exclude positions for which 
the extinction is not defined in the high A v regions, con- 
tain no CO data or have statistical errors for the 13 CO col- 
umn density in excess of 5x 10 16 cm -2 . Such regions account 
for 5.4% and 1.6% of the respective areas in Orion A and 
Orion B. For a given position on the sky, r, the total hydro- 
gen column density, Nh{v) = Nhi(f) + 2Nh2{v), is derived 
from the infrared-derived exti nction using the conve rsion, 



A^(r) = 1.88xl0 21 A,(r) cm" 2 ijRachford et al. 1 120091. The 



total mass of the surveyed area is 



dnN H (r 



(9) 



where mn is the mass of the hydogen atom, fi is the mean 
atomic weight to account for helium (1.36), D is the distance 
to the cloud (410 pc for Orion), and dQ, corresponds to the 
solid angle of an image pixel. The atomic gas column density 
for each position is determined from the velocity integrated 
Hi 21 cm line brightness temperature, T2i(r, v), 



N H i(r) = 1.82 x 10 18 



2i(r,v. 



-Nai.ref cm 2 (10) 



where Nni,ref = 1-1 x 10 21 cm" 2 is the atomic gas column 
density derived from the Hi 21 cm line within the reference 
fields used to establish the mean colours of background stars 
in the calculation of A v . Therefore, the column densities and 
masses derived from A v and the Hi 21 cm line emission are 
relative to column density values in these reference fields. 
The Hi 21 cm line emission is integrated over the velocity 
interval -5 to 21 km s" 1 to account for the more extended 
velocity dispersion of the atomic gas relative to the molecu- 
lar component. The equivalent extinction of this Hi column 
density is A v .hi = iVjjj/1.88 x 10 21 magnitudes. The corre- 
sponding mass of atomic gas, including helium, is 

Mm = /imnD 2 / dQN H i(r) (11) 

J cloud 

For the molecular gas, two components of the mass are de- 
rived from the visual extinction. The first component is de- 
termined within the low A v regime by considering only those 
lines of sight with extinctions in the range, A v .hi < A v < 
A V: 3, corresponding to where 13 CO is underabundant. The 
molecular hydrogen mass within this regime, Mh2,u, is 



Mr 



xlO 21 / dQAv{r)- / dQN HI (r) 



(12) 



for A v ,hi < A v < A Vt 3. The second mass component, 

Mh2,SS, 



M H 2,ss = Lim H D 



1.88xl0 21 / dClAv(T)- / dfiiVHj(r) 



(13) 



for A v > A v> z, corresponds to the mass within the strongly 
self-shielded layers of the cloud. Owing to the large differ- 
ence between the beam size of the HHl 21 cm line observa- 
tions (36') and the pixel size of the extinction maps (1.8'), 
we assume a smooth distribution of the atomic gas column 
density between these angular scales. 

The molecular hydrogen mass is also estimated from the 
distribution of 13 CO column density. Traditionally, one as- 
sumes a constant value of the 13 CO to H2 abundance ra- 
tio throughout the cloud. Here, we assess the errors in the 
13 CO derived mass, M lte using the s tandard 13 CO to 
H2 abundance 2xl0~ 6 (Dic kmanl 1 19781) that is typically 
applied in studies of G MCs (|Carpenter. Snell. fc Schloerbl 
1995; iHever et al. 1 12009) and compare the result to the A v 
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derived values in both the underabundant and strongly self- 
shielded regimes. 

M LTE ,u = M"iH2£> 2 (5 x 10 s ) J dQN( ri CO){r) (14) 
for A v mi < A v < A v * and 



(imniD 2 (5x 10 5 ) / dfWV( 3 CO)(r) (15) 



for A v > A v ,3. Finally, we calculate the 12 CO luminosity in 
each cloud and partition for each A v regime, 



Lco,u = D J dtt J dvT 12 , m b(r,v) 
for A v hi < A v < A v ,3 and 



^co,ss 



dVt \ dvT 12 , m b(r,v) 



(16) 



(17) 



for A v > A v .3. ri2, m t,(r,w) is the 12 CO spectrum for posi- 
tion, r within the map. The derived masses for the Orion A 
and B clouds and the Orion A partitions are listed in Table 2. 
The corresponding CO luminosity to molecular mass conver- 
sion factor, M/Lco = (Mh2,u+M H 2,ss)/(Lco,u+Lco,ss) 
is also calculated. 

The atomic masses listed in Table 2 are limited to the area 
covered by the CO observations and derived from atomic gas 
column densities from which NHi,ref has been subtracted. 
The residual masses reflect the amount of atomic material 
that resides within the foreground and background layers 
of the Orion molecular cloud as well as the small compo- 
nent of atomic hydrogen expected to be present within the 
deep, strong ly self-shielded interior of the cloud owing to 
cosmic rays ((Goldsmith fc Lill2005h . For the lines of sight 
in the vicinity of the molecular cloud over the limited area 
shown in Figures 1 and 2, this residual atomic gas provides 
a small contribution (2-10%) to the total hydrogen column 
density. However, if one were to restore the layer of atomic 
gas represented by Nni,ref, the mass of atomic gas in the 
observed CO area increases to 1.7xl0 4 Mq for Orion A and 
and 2.0 x 10 4 Mq for Orion B. Moreover, the atomic envelope 
of the Orion cloud extends well beyond the surveyed area so 
the total atomic mass of the Orion region is likely much 



larger than these calculated values. iGrenier et al. I (|2005l ) 



tabulated the gas masses in the Orion cloud over a much 
larger area and found a ratio of 12 CO derived molecular to 
atomic gas mass of 1.2:1 that indicates an extended compo- 
nent of atomic gas well beyond the surveyed areas of Figures 
1 and 2. 

Inspection of Table 2 shows that the 13 CO LTE column 
density with a constant 13 CO to H2 abundance of 2xl0~ 6 
provides a reasonably accurate measure of the H2 mass re- 
siding within the self-shielded interiors of the Orion cloud 
despite the large abundance variations within the cloud. In 
this regime, the 13 CO mass slightly overestimates the mass 
determined by extinction by 25% in Orion B. In Orion A, 
the respective masses are identical. The root mean square 
deviation is 23% for the 9 paritions in Orion A with a max- 
imum deviation of 40% occuring within Partition 6, which 
also exhibits an enhanced 13 CO to H2 abundance ratio. The 
self-shielded layer of the cloud is critical to the star forma- 
tion process as it is the substrate from which dense filaments 
and cores develop and from which newborn stars emerge. 



However, this self-shielded zone comprises a fraction of the 
total molecular mass of the cloud - 60% in Orion A and 47% 
in Orion B. 

The 13 CO derived masses with a constant abundance sys- 
tematically underestimate the H2 mass within the low A v , 
13 CO underabundant regime by factors as large as 3. Fig- 
ure !14l shows a cumulative histogram of the inferred H2 mass, 
respectively traced by the infrared-derived A v and 13 CO, 
which resides within varying zones of visual extinction For 
each cloud approximately 50% of the H2 mass resides within 
lines of sight with A v < A„,3 in which the 13 CO abundance is 
low and changing rapidly. While the small opacity of the low 
J transitions of 13 CO are appealing to discer n cloud struc- 
ture and kinematics l|Goldsmith et al~l 120081 ; iHever et al. I 
2009), investigators must recognize that these emission lines 
are insensitive to the more extended reservoir of molecular 
gas that may comprise a significant fraction of the cloud 
mass. In regions of high UV radiation such as partitions 7-9 
in Orion A, this missing fraction can be even larger. 
Dark gas describes the low A v regime in the vicinity of 
molecular clouds in which molecular hydrogen is abundant 
due to efficient self-shielding but the abundance o f 12 CO is 
low (|Grenier et al. 1120051 ; IWolfire fc McKee II2OI0I ). In these 
conditions, 12 CO is underluminous with respect to its stan- 
dard relationship with H2 at moderate values of A v . In the 
A v range to which our maps are sensitive, there is little ev- 
idence for dark gas throughout most of the Orion A and 
Orion B clouds. The M/Lco ratios in the underluminous 
regime are comparable to those derived in the strongly, self- 
shielded regions. Dark gas is more evident in Orion A parti- 
tions 7-9 where the M/Lco ratios are larger than the values 
in the self-shielded regime. In these regions, the UV field is 
enhanced by the UV emitting young stars of the OB associ- 
ation as illustrated in Figure [TT1 allowing for the destruction 
of CO molecules while maintaining a large fraction of molec- 
ular hydrogen gas. 

The conversion factor, Xco, relates the integrated 12 CO 
emission, Wco, to the column density of H2 molecules, 
Ah2 = XcoWco- The use of the optically thick 
12 CO line as a tracer of column density has been 
j ustified by both simple model s of interstellar clouds 
l|Dickman. Snell. fc Schloerbl |l986) and more recently, by 
computational simulations that include molecular c hemistry 
llGlover fc Mac Low Il201ll; iNaravanan et al. II2012TI. 
iDigel. Hunter, fc Mukheriee I l|l995l ) and lDigel et~al~l (Il999f ) 
studied the Xco factor specifically within the Orion molec- 
ular cloud using MOO MeV 7-rays, 12 CO and Hi 21 cm 
line emission spanning a much larger area than the more 
limited imaging shown in Figures 1 and 2. They found 
A"co = l-35x 10 20 cm~ 2 /(K km s _1 ) corresponding to a 
M/Lco ratio of 2.9 M©/(K km s" 1 pc 2 ). This conversion 
includes contributions from helium atoms to the mass bud- 
get. The Ori on Xco value is lower th an the Galactic mean 
using 7-rays l|Strong fc Mattox"lll996r i and much lower than 
a recent determination using dust emiss ion from clouds in 
the solar neighborhood (|Ade et al. Il201l1 ). This gamma-ray 
derived value of Xco is similar to those calculated for Orion 
A (3.2) and Orion B (2.9) in Table 2 using our 12 CO mea- 
surements and A v . 

More recently, the Orion cloud complex was observed b y 
the Fermi Large Area Telescope (|Ackermann et al~ll2012l ). 
This study found a difference in Xco values between the 



© 0000 RAS, MNRAS 000, 000-000 



CO Abundance Variations in the Orion Molecular Cloud 19 



< 

v 

<r 



1.0 



0.5 



0.0 



Orion A 




10 20 
A v (mag.) 



30 



10 20 30 

A v (mag.) 



Figure 14. A cumulative histogram of molecular mass in Orion A and Orion B for A v < A v ' as traced by A v (solid line), 12 CO, 
(dashed-dotted line) assuming a CO luminosity-to-mass factor tabulated for each cloud in Table 2, and N( 13 CO) (dashed line) with an 
13 CO-to-H 2 abundance of 2.0xl0" 6 . 50% of the H 2 masses reside within the low A v regimes of the clouds. 



high longitude end of the Orion A cloud (I > 212) where 
X C o = 2.3xl0 20 cm- 2 /(K km s" 1 ) (M/Lco=4.9) and 
the remainder of the Orion A cloud as well as Orion B 
where X CO = 1.3xl0 20 Cm - 2 /(K km s" 1 ) (M/L 00 =2.8). 
They attribute these differences to the rapid changing CO 
abundance in the diffuse gas regime at high longitudes of 
the Orion cloud and the larger contribution by "dark" H2 
material to the total molecular mass. While our observa- 
tions do not extend as far as the gamma ray imagery, 
we derive an elevated mass weighted average M/Lco of 
3.5 M0/(Kkms _1 pc 2 ) within partitions 1 and 2 for which 
I > 212 compared to the remainder of the Orion A cloud 
where M/Lco=3.1 M /(Kkms~ 1 pc 2 ). We also note the 
high values of M/Lco in partitions 7,8,9 that are close to 
the central region of the OB associaton and illuminated by 
local UV radiation fields. 

Our estimates of Xco in Orion are sm aller than the 
value (2.8xl0 20 cm _2 /(K kms" 1 )) derived bv lParadis et al. 
(|2012T ) that is also based on infrared extinctions using 
2MASS. iParadis et all |2012T ) examined a more extended 
area than that sh own in Figures [I] and [Hand employed a dif- 
ferent algorithm ( Doba shil 201lf ) to account for foreground 
star contamination, which may account for the measured 
differences. 

Our results for the Orion molecular cloud and those de- 
rived from gamma ray measurements emphasize that there 
are large variations of Xco in the interstellar medium. The 
application of the Galactic averaged value to an individual 
cloud can lead to 100% errors in the molecular mass. More- 
over, the scatter in the Wco~A v relationship (see Figures [TCT1 
and Hip demonstrate that the varying environments within 
a molecular cloud can lead to even larger errors in the H2 
column density. 



6 CONCLUSIONS 

In this study, we derive two independent views of the molec- 
ular mass distribution in the Orion A and Orion B molecu- 



lar clouds using 12 CO and "CO line emission and IR stellar 
photometry from 2MASS. These data are analysed to inves- 
tigate the variation of 13 CO abundances over the changing 
physical conditions within the cloud and to evaluate the ef- 
fectiveness of 13 CO emission as tracer of molecular mass in 
interstellar clouds. We find the following results. 

• Three distinct regimes are identified in the N( 13 CO)- 
A v relationship: the photon-dominated envelope at A v < 3 
mag., the fully self-shielded cloud interior at 3 < A v < 10 
mag. in which the 13 CO abundance is stable, and high col- 
umn density zones ( A v > 10 mag.) in which CO can deplete 
onto cold dust grains. The degree to which each regime is 
present depends on the local physical conditions - in partic- 
ular, the strength of the UV radiation field and the temper- 
ature of dust grains. 

• In Orion A, the measured depth of the PDR envelope 
increases with proximity to the centre of the Orion OBI 
association, which provides the primary source of UV radi- 
ation. This spatial dependence emphasizes the requirement 
of a self-shielding layer to develop significant 13 CO abun- 
dance values. The relative 13 CO abundance within the fully, 
self-shielded regime varies across the cloud from a minimum 
of 1.4xl0" 6 to 3.4 xlO" 6 . The maximum 13 CO abundance 
occurs near the Orion Nebula Cluster. 

• The presence of CO depletion onto dust grains depends 
on the local gas temperature that serves as a relative mea- 
sure of the dust temperature. The observed dependence 
points to the role of thermal evaporation of CO molecules 
from grain surfaces whose temperature is above a critical 
value. The regions closer to the Orion OBI association show 
less depletion, which may account for the higher 13 CO abun- 
dances found there. 

• Assuming a constant 13 CO to H2 abundance of 2x 10 -6 
throughout the cloud, the 13 CO, LTE derived column densi- 
ties provide an accurate accounting of the H2 mass in layers 
of the cloud that are fully self-shielded. In Orion, this layer 
comprises ~50% of the total molecular mass of the cloud. 

• Using dust extinction as a measure of H2 column den- 
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sity, we confirm the reduced value of Xco in Orion relative 
to the Galactic mean value previously derived from gamma- 
ray measurements. This difference attests to variations of 
Xco in the interstellar medium. 
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